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Classical properties of low-dimensional conductors: 
Giant capacitance and non-Ohmic potential drop 
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Electrical field arising around an inhomogeneous conductor when an electrical current passes 
through it is not screened, as distinct from 3D conductors, in low-dimensional conductors. As 
a result, the electrical field depends on the global distribution of the conductivity a(x) rather 
than on the local value of it, inhomogeneities of o~(x) produce giant capacitances C(u>) that show 
frequency dependence at relatively low u>, and electrical fields develop in vast regions around the 
inhomogeneities of a(x). A theory of these phenomena is presented for 2D conductors. 
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Quantum confinement of electrons in low-dimensional 
(LD) conductors results in dramatic suppression of the 
electrical screening that changes from the exponential 
in 3D to the Stern power-law screening in 2D Q. It 
manifests itself, e.g., in the softening of plasma modes 
Screening phenomena should also strongly influence 
the passage of currents through the circuits including 
inhomogeneous LD conductors. It is well understood 
now that the applicability criteria of the four- and two- 
terminal Landauer-Biittiker equations || are controlled 
by the electrical screening Q . The latter is also of impor- 
tance for calculating spin-injection transients [H and for 
some different problems of spintronic devices M]. How- 
ever, calculating the effect of screening on the potential 
distribution around a long resistive LD conductor with a 
potential moderately growing at infinity needs the tech- 
niques that differ from the standard techniques used, 
e.g., for plasmons. Meantime, in addition to the hori- 
zontal transport in narrow quantum wells (remarkably, 
few-nanometer thick strips have been created) and in 
single- wall nanotubes electrical transport through 

a chain of individual gold atoms j^] and even through in- 
dividual molecules |nj has been reported. Under these 
conditions developing a regular procedure for calculating 
the effect of screening on the parameters of extended LD 
conductors becomes an important issue. 

In regular 3D conductors the screening length l^o is 
small compared to all geometrical dimensions. There- 
fore, the inhomogeneous electrical field produced in vac- 
uum (or an insulator) around the conductor by the po- 
tential drop inside it is completely screened in the narrow 
surface layer about I^d and does not penetrate into the 
bulk of the conductor. As a result, there exists the usual 
local Ohmic relation between the field E(r) and the cur- 
rent density j(r), E(r) = j(r)/tr(r), where <r(r) is the 
conductivity (the mean free path ^ m .f. p . is assumed to 
be small enough). This is not the case for LD conduc- 



tors with long-range screening. Therefore, the relation 
between j(r) and E(r') becomes nonlocal. 
We show that for LD conductors: 

(i) the potential drop is non-Ohmic because the field 
E(r) is controlled by the global distribution of a(r) rather 
than by its local value, 

(ii) there is no macroscopic limit for the potential drop 
per unit length for a long LD conductor, and 

(iii) inhomogeneous "Ohmic" LD conductors acquire a 
giant capacitance that shows dispersion at low frequen- 
cies and can be measured in ac experiments. 

To make the essence of the phenomena most clear we 
choose the simplest model of a conductor in the diffu- 
sive regime when a local constitutive equation j(r) = 
er(r)VC(r) relating j(r) to the gradient of the electro- 
chemical potential £(r) holds. The 2D conductor resides 
in the sz-plane, its 2D conductivity <r(r) = o~(x), and 
the current j(r)||x. The electrical potential ip(x,y) sat- 
isfies the Laplace equation Atp(x, y) = for y ^ and 
is continuous at the conductor boundary. The boundary 
condition follows from the equation 



n{x) = ep[C(x) + tp(x)], ip(x) = ip(x, y = 0), 



(1) 



that is fulfilled inside the LD conductor in the linear ap- 
proximation in the current j{x). Eq. ([!]) relates (,{x) and 
<p(x) to the nonequilibrium 2D electron concentration in 
the conductor, n(x), and the 2D density-of-states in it, 
p, that we assume to be x-independent. For the upper 
half-plane, y > 0, the boundary condition for the Laplace 
problem, according to the Gauss theorem, reads as 



l<Py(x) ~ <p(x) = C{x), <p' y (x) = d y ip(x, y)\ y . 



(2) 



where I = e/2ire 2 p is the 2D screening length, e being 
the dielectric constant of the insulator surrounding the 
2D conductor. The ratio I/I^d ~ V r B kp/kpw, where tb 
is the Bohr radius, kp is the Fermi momentum, and w 
is the thickness of the 2D layer (in the y-direction) . In 
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a quantum conductor this ratio is large, I/I3D 3> 1, be- 
cause rskp > 1 for a metallic conductor, and Icfw <C 1 
under the conditions of a strong 2D confinement. 

Therefore, the Laplace equation should be solved with 
the mixed boundary condition of Eq. (|J). Let us intro- 
duce an auxiliary function 

y r°° at) dt 1 r°° at) 

(3) 

z = x + iy, that is harmonic for y > and satisfies the 
Dirichlet boundary condition ip{x, 0) = ((x) provided the 
integral exists. Such an approach is advantageous com- 
pared to the standard Fourier transformation technique 
because it works for functions £(x) moderately increasing 
at infinity. For £(2) oc |x| Q sign{x} the integral in Eq. (||) 
converges for a < 1 and diverges as a power of the lead 
length for a > 1. The case of Ohmic leads, a — 1, is 
marginal and results in the logarithmic singularities dis- 
cussed below. Then the function 

/°° dt 
r°° rli 

= -jf ^(x,y + t)e-^j (4) 

is also harmonic in the upper half-plane, y > 0, and sat- 
isfies the mixed boundary condition of Eq. (||) . Because 
C( x ) = 3 l°~{ x ) f° r an one-dimensional flow, the function 
C(x) is known and Eq. (6) provides an explicit solution 
for the potential (f(x,y). 

To find out the basic regularities of the potential dis- 
tribution inside LD conductors (and around them) let us 
consider several model systems. For a resistive contact at 
x = in a perfect conductor, l/a(x) = <5(x)/£, the elec- 
trochemical potential equals ((x) — (j/£)0(:r), where 
S _1 is the contact resistance and Q(x) is the Heaviside 
function. It follows from Eqs. (||) and (0) that 

fp(x,y) = (j'/ttS) arccot(-x/y), 

<p{x) = -^e^ + r 2 e '] dT . (5) 

As distinct from a x )> the potential tp(x) is continu- 
ous near x = but acquires a singular contribution 
(jx/7rE)ln(|x|/Z) 0. Therefore, the electrical field di- 
verges near x = as E x {x) ~ (j /ttIH) ln(|x|/i) inside the 
conductor and in the vicinity of it. For \x\ >> I the po- 
tential ip(x) approaches its asymptotical values slowly, as 
fix) « — Q{x) + jl/irHx. The total drop in ip(x) is equal 
to the drop in C( x )- However, while the drop in C( x ) 
occurs abruptly at x = 0, the potential tp(x) changes 
gradually at the scale of I. 

With the ip(x) found above, the current induced 
nonequilibrium electron concentration n(x), see Ec;. (Q), 



behaves as n{x) oc 1/x for \x\ ^> I, hence, the dipole 
moment of this charge distribution diverges. The giant 
dipole moments developing near the resistive elements of 
the circuit, and electrical fields around them, are an ex- 
ceptional property of LD conductors. In 3D the screening 
is exponential and develops at the scale of I^d- 

If a resistive conductor of the length L and the con- 
ductivity cr(x) = const is connected to two perfect con- 
ductors, then a x < 0) = 0, ((x) = jx/a for < x < L, 
and Ci x ) = jL/<r for x > L. Performing integrations 
in Eqs. (j|) and (^) results in two scales in tp(x), I and 
L. When L 3> I, the behavior of <p(x) at the scale L is 
controlled by the term — (jl/ira) ln(|L — x|/|x|). It has 
the magnitude of about (jl/wa) \n(L/l) and describes the 
electrical field E x (x) that penetrates into the perfect con- 
ductor (and the insulator around it) as deep as by several 
L. Fast changes in (f(x) at the scale / are of importance 
in the vicinity of the contacts, i.e., near x — and x = L. 
However, because the contacts are supposed be perfect, 
both tp(x) and E x {x) are continuous near them. 

In both model systems discussed above, the Coulomb 
nonlocality characteristic of LD conductors manifests it- 
self quite distinctly. First, j(x) and E x (x) are correlated 
at the scales of I and L that may be much larger than 
the scales i m .f.p. and I3D typical of 3D conductors. Sec- 
ond, because of the logarithmic terms the problem of the 
potential distribution in a long homogeneous LD conduc- 
tor has no macroscopic limit. It is only due to the 3D 
contacts fixing the potential drop that the total drops 
in ip and £ are equal, see below. Third, the electrical 
fields developing near the inhomogeneities of the conduc- 
tor deeply penetrate the insulator surrounding it (or vac- 
uum); they can be detected by a scanning microscope. 
These fields also result in the mechanical interaction of 
two LD conductors (their attraction or repulsion) and in 
their mechanical deformation. This mechanism provides 
a macroscopic approach to the current-induced forces in- 
vestigated recently by first-principles calculations fH| . 

The change in E x (x) in response to the modulation of 
o~(x') in a remote point x' in the j= const regime, as 
well as the current induced distributions of E x (x) and 
n(x) around the inhomogeneities, including the field in 
vacuum near the conductor, are the predictions of the 
theory that can be checked experimentally. 

We are now in a position to investigate a realistic sys- 
tem of an inhomogeneous 2D conductor with a conduc- 
tivity a(x) connected by perfect contacts to two classical 
conductors in the planes x — and x — L, their elec- 
trochemical potentials being £(0) = and ((L) = (l, 
respectively. Subtracting xQl/L from C( x )> we arrive at 
the function 

C(x) = C(x) - xa/L (6) 

that equals zero at both ends of the interval (0, L) and, 
therefore, can be expanded in a Fourier sine series 
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Here integration by parts was performed. Substituting 
fT Eq. (fjj) into Eq. (0) and (||) results in the final expres- 
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sions for tp(x,y) and <p{x,y) in the upper half-strip: 
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The potential y) of Eq. (|^) does not depend on y in 
the planes x = and x = L, hence, the potential drop 
in the classical conductors |l3| being small compared to 
the potential drop in the LD conductor is disregarded. 

It follows from Eq. (^) that the typical penetration 
depth of the electrical field into the insulator is about 
L. However, I controls the magnitude of the field that is 
reduced by the factor about L/l when nl/L^> 1. 

Eq. allows to find the kernel K(x,u) = 

8E x (x) / 8C, 1 (x) relating the response of the electrical field 
E x {x) = —<p'(x) at the point x to the variation of the 
resistivity l/a(u) — C'( u )/j a ^ the point u. In 3D this 
relation is nearly local in the diffusion approximation, 
K(x, u) w S(x — u) with the width of about I^d- Differ- 
entiating Eq. (||) with respect to x, taking into account 

that £l = J C,'(x)dx, and performing the variational 
derivative, one finds 



12 cos(nnx / L) cos(Trnu / L) 
K{x,y) = - + -^ 1 + nnl/L ' 

n—l ' 



(10) 



This kernel is completely determined by the geometry of 
the system and does not depend on the specific form of 
a(x). The relation K(x, u) du = 1, providing for the 
equal drop of the electrical and electrochemical poten- 
tials across the specimen, follows from Eq. ( |l0| ) because 
the integral from the second term vanishes. 

Using the identity (1 + nnl/L)^ 1 = J °° eiAexpf— A(l + 
nnl/L)], the summation in Eq. (|T^) can be performed. 
Finally K(x, u) = l/L + f(x — u) + f(x + u) where 
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f(u) = — 



: dx c _ a cos(7tu/L) - exp(-7rlA/L) 
cosh(7rZA/L) — cos(ttu/L) 



Integration in Eq. ( pT| ) over the region A > u/l results 
in a logarithmic singularity in K(x,u) at x w it that is 
much weaker than the S(x — it) singularity in 3D. When 
the LD conductor is long enough, I < L, this singular- 
ity is of the scale K(x,u) ~ (1/0 ln(Z/|a; — it|), and is 
followed by the intermediate power-law asymptotic of 
K(x,u) ~ l/\x — u\ in the region I <§C \x — u\ <C L 
(if / <C L). Under these conditions the overall scale of 



K(x, u) is about l/L, however, for \x — u\ <j I it reaches 
the scale l/l and can exceed considerably the constant 
contribution l/L. In the opposite limit of a short LD 
conductor, L <C I, the shape of K(x,u) depends only on 
L while its magnitude is proportional to l/l. Near the 
singularity K{x,u) ~ (l/V)\n{L/\x — u\). In this limit 
the screening comes mostly from the classical electrodes, 
and the term l/L dominates the kernel K(x,u). 

Electrical field developing in the insulator around in- 
homogeneities of u(x) results in the electrostatic energy, 
and hence in the capacitance of LD conductors Jl4|] . We 
will show that for LD conductors C has a giant mag- 
nitude and shows a frequency dependence at relatively 
low frequencies to. For comparison, in 3D the capaci- 
tance can change at the scale of the geometric capac- 
itance C g — e/AitL (per unit area) when the modula- 
tion of ct{x) is deep enough. However, a considerable 
frequency dependence of C is expected only for frequen- 
cies larger than the inverse Maxwellian relaxation time, 
k> ^ ^>3D = ^ctsd/s, where a^d is the 3D conductivity. 
With uzd ~ cr/w, one finds uj^d ~ ^Tra/ew. Measuring 
the low-frequency dependence of the capacity of a LD 
conductor might become a practical tool for detecting 
inhomogcncities of a(x). 

For an external potential depending on the time t 
as exp(— iuit), the continuity equation reads j'(x) + 
iuien{x) = 0. Eq. ([!]) with £(x) found from Eqs. (||) 
and (Q) and ip(x) found from Eq. (^|) should be used for 
n(x). Finally, the continuity equation reads Eq] 
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d[a(x)C(x)}/dx 
cue sin(7rnx/L) 



ttL <H 1 

n—l 



iml/L 



C(u) cos(nnu/L) du = 0. (12) 



This is an intcgro-diffcrcntial equation for ('(x) that for 
replaces the equation C( x ) — j/ a ( x ) valid in the 
dc regime. The capacitance C{uS) and the active resis- 
tance R(lu), per unit width in the z direction, that can 
be found from Eq. ( |l2| ) depend on the frequency u>. The 
characteristic frequency ujid ~ (J /el is controlled by the 
2D screening length I, hence, UJ2D *C u>3d- The dispersion 
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of C(u) and R(u>) for uj > ti^D should allow to distinguish 
the capacity of a LD conductor from the geometrical ca- 
pacitance. 

The capacitance C(uS) is related to the total current in 
the circuit J(t) including both j(x,t) and the displace- 
ment current jd{x,y,t) originating from the charges ac- 
cumulating at the electrodes. In a quasistationary regime 
the current J(x, t) integrated over y is independent of x 

J{t) = o(x)dQ(x, t)/dx + jd(x, t) — const. (13) 

The current jd{x, t) can be found from the second term of 
Eq. (^) that represents the contribution to f{x, y) coming 
from the area adjacent to the LD conductor: 



j d {x) = {ius/An) / dy dip(x, y)/dx 

J —OC 

= -^g i tl (l+J/L) j D ^C»cos( — ). (14) 

In the lowest order in w, one can apply the static equa- 
tion C( x ) = Cl/ Ro&{x) and express jd{x) in terms of the 
cosine- Fourier components of the resistivity a^ 1 (x) 



Rn = I dx cos(irnx/L)/a(x), n>0. (15) 
'o 



Therefore 



iuje Cl Rn cos{imx/L) 

n—1 



(16) 



Substituting Eq. ( [l~6| ) into Eq. (|i~3|), dividing it by a (a 
and integrating over x results in 



JRo = U 2~5~ ~7T~ 

7T Ro , nil 



Ro Z - 1 nil + nrd/LY 



(17) 



Therefore, the final expression for the low-frequency ca- 
pacitance (per unit width in the z direction) is 



C = 



-Y — 

7T 2 <H nil 



R n /Rl 



' n(l + iml/L) 



(18) 



Eq. jlq) reveals a number of remarkable properties of 
the capacitance. First of all, C depends on the conductiv- 
ity <j(x) only through the ratios (R n /Ro) 2 [i.e., through 
the degree of the inhomogeneity of the LD conductor], 
is independent of the total resistance Ro, and vanishes 
for a homogeneous conductor Jig . Second, C — > when 
I — ► oo, hence, it comes exclusively from the screening 
in the LD conductor. Third, if <j(x) is scaled as a(x/L), 
the ratios R n / R$ do not depend on L while C increases 
monotonically with L and saturates as L — > oo. 

The most remarkable property of C is its giant mag- 
nitude under the conditions when I < L and the inho- 
mogeneity is strong enough, J^nLi R n / n ~ -^o- Indeed, 



in this case the ratio C/C g ~ L. This implies that a 
2D conductor with the thickness w acquires an effective 
thickness L ^> w. This increase originates from the pen- 
etration of the electrical field into the insulator at the 
depth of about L. Only in the opposite limit, I ^> L, 
the ratio drops to C/C g ~ L? jl because the field in the 
insulator is reduced by the factor L jl as it is has already 
been mentioned above. 

Giant capacitances are a general property of LD con- 
ductors that manifests itself quite dramatically for a re- 
sistive contact in a perfect conductor discussed above. 
With ip(x,y) of Eq. (||), the x component of the elec- 
trical field in the x = plane equals E x (0,y) — 
— (j/nYtl) exp(y/Z)Ei(— y/l), where Ei(— y/l) is the in- 
tegral exponent function. For y ^> I the field equals 
£^(0, y) ~ j /iryYi. Therefore, the integral over y that en- 
ters into the displacement current jd(x = 0), see Eq. (|lj), 
diverge logarithmically for large y, hence, the capaci- 
tance C that is proportional to that integral, also di- 
verges. This integral is cut-off by the geometry of the 3D 
conductors surrounding such a contact and screening it. 
Finally 



C« (e/27r 2 )ln(L/0 



(19) 



where L is the cut-off length. 

Our approach to the Coulomb problem as developed 
for the diffusive transport is not restricted to it. In the 
classical regime, quantum wires should show properties 
similar to those of 2D systems but the mathematical tech- 
niques are different. As applied to the spin injection, the 
dimensional effects do not influence the two-terminal re- 
sistance controlled by ((x) but should be of importance 
for optics and transients because the Poisson equation is 
involved. 

In conclusion, we have developed a technique for solv- 
ing the Coulomb problem for 2D conductors and have 
investigated their properties including the strongly non- 
local relation between the resistivity and electric field, gi- 
ant capacitances, and electrical fields developing around 
inhomogeneities. 
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